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ABSTRACT. We consider second order differential operators P with polynomial 
coefficients that preserve the vector space 14 of polynomials of degrees not greater 
then k. We assume that the metric associated with the symbol of P is flat and that 
the operator P is potential. In the case of two independent variables we obtain some 
classihcation results and find polynomial forms for the elliptic A2 and G2 Calogero- 
Moser Hamiltonians and for the elliptic Inosemtsev model. 

1. Introduction 


In the paper [Tj a transformation that brings the elliptic Calogero-Moser operator 

3 

H = A + g^p{xi-Xj), ( 1 . 1 ) 

i>j 

where g is arbitrary constant, to a differential operator with polynomial coefficients. This 
problem is easily reduced to two-dimensional case. Indeed, in variables yi = Xi — where 
Xi the Laplacian is given by 


A = -3 


2/^2 52 52 

3 V5y2 Qy'i dyidy2 


The variable Y is separated and we reduce fll.ip to 


« = -3 


52 92 

+ 


92 


dyl dyl dyidy^ 


+ 9^ p{yi - Vj), 


i>j 


where 7/3 = -yi 


X = 


y 2 . It turns out that the change of variables 

p'{yi) - p'{y2) p{yi) - p{y2) 


y = 


( 1 . 2 ) 


(1.3) 


p{.yi)p'{.y2) - p{.y2)p'{.yi) p{.yi)p'{.y2) - p{.y2)p'{.yi) 

and a gauge transformation of the form L —> fLf~^, where / is an appropriate function in x, y, 
bring fll. 2 p to a differential operator of the form 


^2 ^2 52 a 9 

L = a{x,y)— + 2b{x,y)—— + c{x,y)— + d{x,y)— + e{x,y)— + f{x,y) (1.4) 
ox^ oxoy oy^ ox oy 

with polynomial coefficients. We call L polynomial form of fll. 2 p . The set of coefficients a,b, 

and c is called the symbol of L. 

A polynomial form for the elliptic case have been found in [1] by a deformation of known 
[ 2 ] polynomial form of the trigonometric model. Concerning polynomial forms for the elliptic 
Calogero-Moser see |3]. 
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The goal of this paper is a description and an investigation of a class of operators of the form 
fll.dp that contains polynomial forms o many operators related to simple Lie algebras |1]. This 
class is selected by postulating of some properties of polynomial form for operator fll.2p . 

It is clear that for this polynonial form L 

1 : the contravariant metric g defining by the symbol of L is flat □ 
and 

2 : L can be reduced to a self-adjoint operator by a gauge transformation L —> fLf~^. □ 

Furthermore, it can be shown that for some special values of the parameter g in fll.2p 

3 : the operator L preserves the vector space 14 of all polynomials P{x, y) such that deg P ^ n 
for some n ^ 2. □ 

and that 

4: L is invariant with respect to the involution x = x, y = —y. □ 

In Section 2 we consider diffrential operators of second order with one independent variable. 
In this case the assumptions 1-2 are automatically satisfied. We describe all operators that 
satisfy Property 3 . As a result we arrive at the polynomial forms for the Lame operator and 
for the more general Darboux operator [S]. 

Our main goal is to investigate two-dimensional polynomial operators fll.ip that satisfy 
Properties 1 - 3 . In Section 3.1 we describe operators that satisfy Property 3 , define an action 
of the group SL^ on the vector space of such operators, and find some of invariants for this 
action. 

In Section 3.2 we find all operators with elliptic and trigonometric symbols (see Definition 1) 
that satisfy Properties 1 - 4 . As a result, we arrive at polynomials forms for the elliptic A 2 and 
G 2 -Hamiltonians and for the Inozemtsev elliptic model. The A 2 and G 2 polynomial forms are 
equivalent to ones presented in []. The polynomial form for the Inozemtsev model is probably 
new. Possibly it can be deduce from the rational form found in [B]. It would be interesting to 
find the Hamiltonians in the Schrodinger form corresponding to five trigonometric sysmbols 
obtained during the classification. 

By miracle, all known operators L satisfying Properties 1-3 turn out to be integrable (i.e. 
possess differential operators that commute with L). 

It is clear that if a differential operator obeys Property 3, we can find several eigenvalues 
and corresponding polynomial eigenvectors by means of linear algebra. In [7] such operators 
are called quasi-exact solvable. 
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For the elliptic models the flat metric g depends on the elliptic parameter. One of the reasons 
why such a metric could be interesting in itself is that families of contravariant metrics with 
linear dependence on a parameter are closely related to the Frobenius manifolds j8] 

Acknowledgments. The author is grateful to M. Matushko, A. Turbiner, and E. Ferapontov 
for useful discussions. He is also grateful to MPIM(Bonn) for hospitality and financial support. 

2. Differential operators with one independent variable 

Suppose that an order m polynomial operator of the form 

Q = 5^a,(a;)— (2.1) 

1=0 

preserves the vector space 14 of all polynomials of degree 4 n, where n ^ m. 

Conjecture (A. Turbiner [^). Any such operator Q is a (non-commutative) polynomial of 
order m in first order generators 

Ji = 1, J2 = J 3 = J4 = - nx. □ 

ax ax ax 

Notice that the Lie algebra generated by Ji, ..., J 4 is isomorphic to gl{2). 

For small m this statement can be straightforwardly verified using the following 
Lemma 1. If operator (12.1 p preserves the vector space 14, where n ^ m, then the degree of 
the polynomial ai is not greater then m + i. 

Proof. Define a weighting scheme in the ring of differential operators with polynomial 
coefficients by setting weight(x) = 1 and weight ^ = —1. Let Q = Yl^i=-m Qi decomposition 

of Q into the sum of homogeneous differential operators Qi, where weight(Qj) = i. We have to 
show that p ^ m. Indeed, if p > m then Qp[x^) = Qp[x'^~^) = ■ ■ ■ = Qp{x'^~'^) = 0. Since the 
kernel of the differential operator Qp of order m contains m + 1 linearly independent functions, 
we get Qp = 0. □ 

It follows from Lemma 1 that any operator P of second order that preserves 14 has the 
following structure: 

(P d 

P = (a^x^ + asx^ + 02 ^^ + aix + ao)^-^ + {hx^ + b2x‘^ + bix + bo) — + C2X^ + Cix + Cq. (2.2) 

dx^ dx 

Taking into account that P{x^) and P{x^~^) have to be polynomials of degree not greater then 
n, we find that the constants in the coefficients of (12.21) are connected by the following identies 


63 = 2(1 - n ) 04, 


C2 = n{n — 1) 04, 


Cl = 77.(03 - na^ - 62). 


(2.3) 




4 


The transformation gronp 

X —>•-——, P {s^x + S4) '^P{s^x + 54)*^, ( 2 . 4 ) 

S^X + S4 

acts on the nine-dimensional vector space of snch operators. The coefficient a{x) at the second 
derivative is a fonrth order polynomial, which transforms as follows 

a(a;)—>■ (ssx-f 54)^0( 2 - 5 ) 

\ C 3 X * 1 “ C 4 / 

If a{x) has fonr distinct roots on the Riemann sphere, we call the operator P elliptic. In the 
elliptic case nsing transformations (12.Sp . we may rednce a to 

a{x) = 4 a;(x — l)(a; — k), 

where k 7 ^ 0,1 is the elliptic parameter. 

Define parameters rj-i,..., ns by identities 

60 = 2(1 -|- 2ni'), bi = —4^(n -|- l)(ni -|- 1) -|- nn 2 -|- , 62 = —2 (3 -|- 2ni -|- 2n2 -|- 2n3), 

n = + n 2 + n 3 + n^), n 5 = cq + n 2 (l - n 2 ) - ns) (ni ns)^ n(ni n 2 )^. 

Then the operator H = hPh~^, where 

np 712 713 

h = x~{x — l)~{x — k) 2 


has the form 

a'{x) d ni(l-ni)n n 2 (l - n 2 )(l - n) , n 3 (l - n 3 )K(K - 1) 

a[x)— + ^—— + n^ + nA[l-ni)x^ -^-^-^-. 

dx^ 2 dx x X — 1 X — K 

After the transformation x = f{y), where f'^ = 4/(/ — 1)(/ — k) we arrive at 

ni(l-ni)n n 2 (l - n 2 )(l - n) n 3 (l - n 3 )n(n - 1) 

H = ^ + n5 + n4(l - n4) / H- 7 -^- - —^-^- - -. 

dy^ f / - 1 f -K 

Here n* are arbitrary parameters. When n = — ^(ni -|- n 2 -|- ns -|- n 4 ) is a natnral nnmber, the 
initial polynomial operator f|2.2p preserves the finite-dimensional polynomial vector space Vn- 
Another form of this Hamiltonian (np to an additive constant) is given by 


H = ^-Fn4(l-n4) p(a;)-f-ni(l-ni) p(a;-l-a;i)-l-n2(l-n2) p{x+u;2) + n3{l-n3) p{x+uJi+U2), 

where Ut are half-periods of the Weierstrass fnnction p{x). If ni = n 2 = ns = 0 we get the 
Lame operator. For generic case H defines the so called Darbonx-Treibich-Verdier model |10j . 
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3. Two-dimensional case 

3.1. Structure of operators having Property 3. The following statement describes the 
vector space of differential oparators we are dealing with. 

Theorem 1. An operator L of the form fll.4p with polynomial coefficients satisfies Property 
3 iff the coefficients of L have the following structure 

a = qix^ + q2X^y + q^x'^y'^ + kix^ + k2x‘^y + k^xy"^ + aix^ + a2xy + a^y'^ + a4X + a^y + a^; 

b = qix^y + q2x'^y'^ + qsxy^ + | (^k^x^ + {ki + k^)x^y + {k2 + kf^)xy‘^ + k^y'^^ + 

hix'^ + b 2 xy + 63 ?/^ + b^x + b^y + b^] 


c = qix^y‘^ + q 2 xy‘^ + q^y^ + kix‘^y + k^xy"^ + k^y^ + Cix^ + C 2 xy + + C4X + c^y + c^; 

d = (1 - n) (^ 2 (qix^ + q2X^y + qsxy^) + kyx'^ + (/c2 + kg - kQ)xy + k^y'^'^ + dix + d2y + dg] 
e = (1 - n) + q^y^) + A;4X^ + {k^ + h - ki)xy + + eix + 62?/ + 63; 

/ = n(n - 1) (qix'^ + q 2 xy + qgy'^ + {kj - ki)x + {kg - kQ)y'^ + /i. 

(3.1) 

Proof. Define a weighting scheme in the ring of differential operators with polynomial 
coefficients by setting weight(x) =weight(y) = 1 and weight^ = weight^ = —1. Let L = 
'YlT=- 2 ^i decomposition of L into the sum of homogeneous differential operators Lj, 

where weight(Lj) = i. Let us prove that m ^ 2. In other words we are going to show that the 
degrees of polynomials a, b, c in fll.4p are not greater then 4, the degrees of d, e are not greater 
then 3 and the degree of / is not greater then 2. Let 


r 

Lm, — <x ——+ 2/3 


dx'^ 


+ 7: 




dxdy dy 


d d , 

+ ^ Q-^ Q-^ 

OX oy 


If m > 2 then it follows from Property 3 that Lp{x^y^) = 0 for i + j = n,n — l,n — 2. Identities 
Lp{x'^) = Lp{x'^~^) = Lp{x'^~‘^) = 0 imply a = = A = 0. It follows from Lp{y'^) = Lp{y^~^) = 

Lp{y'^~‘^) = 0 that 7 = /i = 0. Now the identity Lp{yx'^~^) = 0 leads to /3 = 0. Thus, m ^ 2. 
It is easy to verify that the conditions L 2 {x^y^) = 0 for i + j = n,n — 1 and Li{x^y^) = 0 for 
i + j = n are equivalent to fl3.ip . The operator L 2 contains 3 parameters in Li we have 8 
parameters ki and Lq, L^i, L -2 are arbitrary. □ 

Remark 1. It is easy to verify that fl3.ip is equivalent to the fact that L is a quadratic 
polynomial in generators 


/ - ^ 
ox 


J - ^ 
dy 


T ^ 

J3 = X — , 

ox 


T ^ 

J 4 = y^, 

OX 


T ^ 

J5 = X—, 

dy 


A) y p, ) 

dy 
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These operators form a Lie algebra isomorphic to s/(3). □ 

Operators Ji in the case of arbitrary number k of independent variables have been presented 
in [ 2 ]. The Lie algebra generated by them is isomorphic to s/(A: + l). Apparently, any polynomial 
operator that preserves the vector space of polynomials Vn for sufficiently large n is a polynomial 
in the generators Jj. Following the line of proofs of Lemma 1 and Theorem 1, one can easily 
prove the following statement toward this conjecture. 

Theorem 2. Suppose differential operator of order m 

L= Y. . 

iiH- 

with polynomial coefficients preserves the vector space of all polyniniials in xi,x*, of degree 
^ k, where k ^ m. Then deg ^ m + A + • • • + 4- D 

The dimension of the space of all operators of the form fl3.ip equals 36. The group SL^ acts 
on this vector space as follows 

i = ^, 9 = L = (3.3) 

where P, Q, R are hrst degree polynomials in x and y. The SL^ representation dehned by 
formula fl3.3p is a sum of irreducible representations with representation spaces Vi, V 2 and V 3 
of dimensions 27, 8 and 1 correspondingly. The vector space Vi is spanned by coefficients of the 
symbol. As a basis in V 2 one can choose 

xi = Sky — k^ — 7ki, X 2 = 5ks — k 2 — Tk^, X 3 = 5di + 2(n — l)(2ai + 62 ), 


X4 — 5ei + 2(n — l)(26i + C2), X5 — 5 d 2 + 2(n — 1)(2&3 + 02), Xg — 5e2 + 2 (ji — l)(2c3 + 62)) 

X7 = 5^3 + 2(n - l)(a4 + 65), Xg = beg + 2 {n — l){b4 + cg). 

The generic orbit of the SLg action on V 2 has dimension 6. Two polynomial invariants of this 
action are given by 

/l = Xg - X3X6 + Xg + 3X4X5 + 3 (? 7 , - l)(XlX 7 + X2X8), 


and 


I2 = 2X3 - dXgXg - 3X3Xg + 2 Xg + 9X4X5 (X3 + Xg) + 


9 (n — l)(XiX3X7 + X2XgX8 — 2XiXgX7 — 2 x 2 XgXg + 3X2X4X7 + dXiXgXg). 
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Assume now that an operator L given by fl3.1l) obeys Property 4. Then the coefficients of L 
have to satisfy the additional symmetry conditions 

a{.x,-y) = a{x,y), b{x,-y) =-b{x,y), c{x,-y) = c{x,y), 
d{x,-y) = d{x,y), e{x,-y) =-e{x,y), f{x,-y) = f{x,y). 

This symmetry is not destroyed by a subgroup 

ax + /3 _ y 


X = 


y = 


a5 — (3^ = 1 . 


(3.4) 


7X + h 7X + h 

of the group fl3.3p isomorphic to SL 2 . Transformations L = ciL + C2 and y = Csy are also 
allowed. 

Transformations fl3.4p act on 15-dimensional vector space of coefficients of polynomials a, b, c. 
The representation spaces for irreducible components of this S'L 2 -representation have dimensions 
5, 3, 3, 3,1. Let us write the coefficients a, b and c in the form 


a — P + Qy'^ 1 


h=\(P'-R)yP-Q'y\ 


c ^ S + i—P" - -rC + T)y^ + -Q"y\ (3,6) 


12 


where degP = 4, degQ = degi? = deg S' = 2, and degT = 0. Then the coefficients of these 
polynomials correspond to irreducible components. Namely, under transformations fl3.4p each 
of these polynomials changes by the rule 

^(x) = + (3.6) 

V 7 X -f 0 / 

Definition 1. Symbol (13.5p of operator L is called elliptic if the polynomial P has four 
different roots on the Riemann sphere. It is called trigonometric if P has the only one double 
root. Sometimes we will call operators with elliptic (trigonometric) symbols elliptic [trigonometric). 
□ 

In the elliptic case we may reduce P{x) to one of the following forms: 


P[x) = x(x — l)(x — K.) or P[x) = (x^ — l)(x^ — k) 


by transformations (13.4p . In the trigonometric case we may assume without loss of generality 
that 

F(x) = x(x — 1) or that P{x) = (x^ — l)x^. 

Example 1. For any constants a 7 ^ 0, /5, A trigonometric operator (II. 4p with coefficients 

a = x^(x^ -I- y'^) -f ax^ -f I3y‘^, b = xy[x^ + y"^) -|- (a — f3)xy, c = y‘^[x^ + y‘^) + (3x‘^ -1- ay"^, 


d = 2[n — l)x(A — x^ — y^) 


e = 2[n — l)y[X — x^ — j/^), / = n[n — l)(x^ -f- y^) 














satisfies Properties 1-4. The operator possesses the discrete group of symmetries isomorphic to 
D 4 , generated by the reflections 

X —)■ —X, y ^ y, X —)■ X, y —)■ —y, x ^ y, y ^ x. (3.7) 


3.2. Classification of fiat symbols. The contravariant metric associated with an operator 
of the form (ll.4p is defined by 


= a, = b, = c. (3.8) 

This metric is flat iff 7?i,2,i,2 = 0, where 7?i,2,i,2 is the only non-trivial component of the curvature 
tensor. This is equivalent to an identity of the form 

21 b'^ttxx — ‘^abbxx + a^Cxx + ^bcaxy — 2 ( 6 ^ -|- ac)bxy + 2abcxy+ 

^ S (3.9) 


(?ayy — 2bcbyy + b'^Cyyj X D + flrst order terms = 0, 

where we denote by D{x,y) the determinant a{x,y)c{x,y) — b{x,yy. In this paper we assume 
that D ^ 0. 

Example 2. The elliptic symbol 

a = (x^ — l)(x^ — k) + (x^ -I- k) y^, b = xy {x"^ + y"^ + 1 — 2k), (3.10) 


c= {k — l)(x^ — 1) -I- (x^ -f- 2 — k) y^ -I- y'^ 

is flat. It defines a linear pencil of polynomial contravariant flat metrics |8]. 

The main observation our classification is based upon is the following 
Lemma 2. Suppose that a,b,c are given by fl3.5p with P(x) 7 ^ 0, 5'(x) 7 ^ 0 and the 
corresponding metric fl3.8p is flat; then any root of the polynomial S' is a root of the polynomial 
P. 

Proof. Substituting fl3.5l) to fl3.9l) and setting y = 0, we get 

P (^ - APS'^P” + 8PSP'S' + S'^P'^+APS'^R’ - APRS S' + SP'^SS" - UP'^S''^ + IbQ^^ -R^S'^^ =0 

It follows from here that S' is a divisor of P‘^S'^ . If the roots of S are distinct this proves the 
statement. If S has a double root, then substituting S = Si, where S'! is a polynomial of first 
order, into above expression, we find that S'! is divisor of P^S'(^. □ 

In the elliptic case without loss of generality we set P(x) = x(x — l)(x — /c). The polynomial 
P has the roots 0,1, cx), k, which can be arbitrarely permuted by some transformations fl3.6p . 
It follows from Lemma 2 that there are two alternatives: A: S has a multiple root, which is 
one of the roots of P; and B: S has two distinct roots. We may put S = ax^ in Case A 
and S' = (Tx(x — 1) in Case B. The constant a can be normalized to 1 by a scaling of y. 
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In both cases seven unknown coefficients of the polynomials Q, R,T can be easily found from 
the system of algebraic equations equivalent to f|3.9p . 

Proposition 1. In Case A the condition = 0 is equivalent to 

5 

P{x) = x{x — l)(a; — k), S{x) = x"^, R{x) = —-(x^ — 2x + 3k — 2kx), 

o 

Q{x) = -(x^ — X + 1 + k"^ — KX — k), T = 0. □ 

9 

Proposition 2. In Case B we get 

P(x) = x(x — l)(x — k), S{x) = x(x — 1), R{x) = —3(x^ — 2kx + k), 

Q{x) = -(x^ — 2kx + 2k^ — k), T = -{2k — 1). □ 

2 3 

It is easy to verify that the symbol from Proposition 2 can be converted into the symbol f|3.10|) 
from Example 2 by a transformation of the form f|3.4p . 

Consider now the trigonometric case P{x) = (x — xo)^(x — xi)(x — X 2 ). According to Lemma 
2 we have the following non-equivalent possibilities: 

a; S'(x) = (x — Xi)^, b; 5'(x) = (x — Xi)(x — X 2 ), 

c: 5'(x) = (x — Xo)(x — Xi), d; S'(x) = (x — Xq)^. 

Without loss of generality we fix xq = 00, Xi = 0, X 2 = 1 or, in other words, P = x(x — 1). 
Proposition 3. In Case a there exists only one trigonometric symbol given by 

P = x(x —1), S = x^, Q = o’ -^=x(3 —2x), T = 0, □ 
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Proposition 4. In Case b we get 


P = x(x — 1), S'= x(x — 1), Q = 1, i? = 3(l —2x), T = —-, □ 

o 

Proposition 5. In Case c we have 

P = x(x-1), S' = X, Q = l:, R = 3, T=^. □ 

Proposition 6. In Case d there exist two families of trigonometric symbols depending on 
arbitrary parameter a. Up to scalings they are given by 

a(x, y) = x(x — 1), b{x, y) = 0, c(x, y) = 1 + ay^ 

and by 


2\2 


a{x,y) = {2x-l) y +a{Ax -Ax-y ), b{x,y) = 2{2x - 1) y {1 + y ), c{x,y) = A{1 + y ) 

□ 
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3.3. Property 2 and linear terms for elliptic symbols. Property 2 means that we can 
reduce L to a potential form 

L = Ag + U 

by a gauge transform L = hLh~^ with a proper function h. Here Ag is the Laplace-Beltrami 
operator corresponding to the metric fl3.8p . It is easy to show that it is possible iff the coefficients 
of fll.4p satisfies the relation 

d / be -cd + c{a^ + by) - b{b^ + Cy) \ 

dy\ D J 

d /bd — ae + a{b^ + Cy) — hi^a^ + by)\ 

dr V D / 

where D = ac — b^ . Properties 1 and 2 guaranty that L can be reduced to the form 

d^ 

by a gauge transform and by a proper (complex) change of variables. 

Notice that given polynomial symbol of operator f|1.4p . fl3.ip . condition fl3.1ip is equivalent 
to a system of algebraic linear equations for the coefficients of polynomials d and e. So, for any 
symbol from Section 3.2 we can easily found admissible d and e. The coefficient / is uniquely 
dehned up to an additive constant by formulas f|3.ip . 

For example, for the elliptic symbol from Proposition 1 condition fl3.1ip leeds to 

d = -{1 — n) ^3(5r^ — Ax — Akx + 3k) + (2r — 1 — k) , 

e = ^(1 — n) y ^9r + y'^ — 6k — 6 ^, f = ^n{n — 1) ^6r + . 

Therefore, in this case we have no arbitrary constants in d and e except for the coupling constant 
n and the elliptic parameter k. 

However there is a non-trivial observation concerning parameters in d and e. Due to the 
involution y —>■ —y the symbol admits the transformation x = x, y = y'^. For polynimial symbol 
thus obtained we can again hnd addmissible d and e. It turns out that now they contain an 
additional parameter decides n and k ! We consider the corresponding model in the next Section 
in details. 


4. A2 and G2 elliptic models 

Applying the transformation x = x,y = y'^ to the symbol from Proposition 1, we get a new 
polynomial symbol with 

a = x{x — l)(a; — k) H—{1 — x + x'^ — k — kx + k^) y, 

9 




















11 


b = ^{7x‘^ — 8x — 8kx + 9k) y + ^(2a; — 1 — r) y'^, (4.1) 

o 9 

4 4 

c = ix'^y H—(4a; — 3 — 3 k) y‘^ + - y^- 
3 9 

The determinant D for this symbol is given by 

D = K, K = ksy^ + Qk2y‘^ + 9kiy + 108A:o, 

where 

ks = {k — 1)^, ko = x^{x — l)(a; — k), 
k2 = {k + l ) x ‘^ + 2(k^ — 4k + l)a; — (k — 2)(k + 1)(2k — 1), 

k^ = x'^ + 8(k + l)x^ - 2(4k2 + 23k + A)x‘^ + 36k(k + l)a; - 27k^. 

The transformation x = x, y = y"^ changes the weighting scheme in the ring of differential 

operators. For the new variables we have weight(a;) = 1, weight(?/) = 2, weight^ = —1 and 

weight^ = —2. Denote by the vector space spanned by monomials x^y^, where i + 2j ^ n. 

Theorem 3. An operator L of the form (11.4p with polynomial coefficients preserves Wn, 
n ^ 4, iff the coefficients of L have the following strnctnre 


a = kix'^ + k2x'^y + fcaa;^ + k^xy + aiy + a2x‘^ + a^x + 04; 

b = 2 kix^y + 2 k 2 xy^ + k^y'^ + k^x^ + k^x'^y + bixy + b 2 X^ + b^y + 64^^ + b^x + feg; 

c = Akix'^y'^ + Ak2y^ + Ak^x^y — A{k^ — kQ)xy‘^ + ciy"^ + C2yx‘^+ 

C'iX^ + Cixy + c^x^ + CQy + c^x'^ + c^x + cg; 

d= {1 — n) {^kiX^ + 2k2xy + k^y^ + kjx"^ + dix + d2] 

e = 2(3 — 2 n) {k\x^y + k2y‘^ — k^xy^ + 2(2 — n) {k^x^ + k^xy^ + 2 k'jxy + eiy + 620;^ + e^x + 64; 


f = {n — l)n {kix^ + k 2 y — k^x^ — n kjx + fi- 


□ 


(4.2) 

We will skip the "spectral"parameter fi in all fnther formnlas. 

Writing L in the form (14.21) . where a,b and c are dehned by (14. ip . we hnd that condition 
(13.111) implies 

1 5 

d = -(1 — n){6x + y){2x — 1 — k) H—(x^ — 2x — 2 kx + 3k), 

9 3 

2 2s 

e = -{9x^ + I2xy + y^ — 9y — 9Ky) + —(3x^ + xy — y — Ky) + 

9 3 

2(n-l),, 2 


9 


(9x — 15xy — 2y + 9y + 9ki/), 
n(n-l) 

/ =-^-(3x + i/) 


9 


s 

- nx. 

3 
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We see that this formulas involve an extra parameter s. 

Introduce parameters m, by identities 

n = —3mi — m 2 , s = 1 + 3mi + 3m2. 

Then 


X 




hLh = Ag + 7712(1 — m 2 ) -h 3mi(l — mi)-h A 


y 

mi m 2 

Here A„ is the Laplace-Beltrami operator, h = K~ , and 


K 


(4.3) 


P = 3x^ — 6(k + 1) x'^ + {y + Ky + 9k) x — 2(k^ — k + l)y, 

A = ^ ^ {3mi + m2)(l + 3mi + 3m2). 

o 

Now we are to hnd a transformation of the form x = (p(i/i, 1 / 2 ), y = 1 / 2 )^ that brings Ag 

to the form (cf. fll.2p l 

1 / 

+ 


A = — 


3 \dyl dyl dyidy2 
Loking for a transformation of the form of formal series 

00 00 

i=0 i=l 

we find several first coefficients of these series in terms of (po and its derivatives, derive a first 
order ODE for (po, and after all write the transformation in the following closed form 


X = 


f{yi)^f{y2) - f{yi)^f{y2) 
f{yi)f'{y2) - f{yi)f'{y2) 


y = -l2 


f{yi)f{y2){f{,yi) - f{y2)) 
f{yi)f'{y2) - f{yi)f'{y2) 


where f'^ = 4/(/ — 1)(k — /). Under this transformation the operator (14.31) becomes 


’^ = -3 


Q2 Q2 

+ 


Q 2 


dyl dyl dyidy 2 


+ V{yi,y 2 ), 


where 


U = (mi - l)mi (^p{yi - 1/2) + p{ 2 yi + 1/2) + p{yi + 21/2)) + 

+ (m2 ^l)m2 _ 

This is just the elliptic G 2 Calogero-Moser model (4]. The elliptic 442-model corresponds to the 
special case m 2 = 0. The invariants of the p-function are related to the parameter k as follows 

92 = - «: + 1), 93 = - 2)(k + 1)(2k - 1). 


The polynomial form of the G 2 -model preserves Wn if n = —3mi — m 2 is a natural number. 
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5 . The Inozemtsev elliptic model 

In this section we investigate the model related to the symbol from Proposition 2 written in 
the form (Id.lOj) . Polynomials d and e for symbol (Id.lOh depend on three arbitrary parameters 
becides n and k. But due to D 4 -symmetry the symbol admits the transformation x = y = 

As a result of this transformation and a scaling we get 

a = x{x — l)(x — k) + (1 — k) x{x + k) y, 
b = x{x + 1 — 2k) y + {1 — k) X y'^, 
c = {1 — x) y + {x + 2 — k) y'^ + {1 — k) y^. 

For this symbol we still may use the anzats of Theorem 1. 

It follows from that 

d = Aix(x + y — Ky) + X2{1 + y — Ky) + p x, 
e = \iy{x + y- Ky) + \^{x -l) + qy, 
f = A4(x + y- Ky) + As, 

where 

Kp + (1 — K)q + Xi( 2 k — 1 ) + A 2(2 — k) A 3 (l — k^) = 0 . 

Thus, d and e depend on hve arbitrary parameters ! Dehne constants m,no,ni, 712,113 from 
relations 


2 A 4 — 3 T 4?ti T 2714 2ti2 T 2713 , 2 A 2 — T 2112 ), 2 A 3 — —1 — 2714 , 

— 2 A 4 = Tio(l + 4m + 2 tio + 2ni + 2n2 + 2ns), 

2p = —1 — 3k — 4:Km + 2 ni — Akui — 2n2 — 2Kn2 — 2ns, 

2q = 4 — 3 k + 4m — 4Kni + 4ni — 2Kni + 2n2 — 4Kn2 + 2ns- 
Let us put L = —4L, where L is the operator (II .41) with the polynomial coefficients dehned 
above. Reducing the operator L to the Schrodinger form by the transformation 

(/(x)-!)(/(?/)-1) 


= f{x)f{y), 


y = 


K — 1 


where 


f^ = 4f{f-l){f-K), 


and by a proper gauge transformation, we get 

3 

H = A + 2 m(m - 1 ) (^p{x + y) + p{x - 2 /)) + X] ni{ni - l){p{x + Ui) + p{y + a;*)). 


i=0 


where a;i,a;2 are the half-periods of the Weierstrass function p(x), Uq = 0 , and ijjs = 0J1 + 0J2- 
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This is so called Inozemtsev BC 2 Hamiltonian HD Its polynomial form preserves I 4 if 

k = ~{2m + 

is a natural number. 


CnUCOK JlHTEPATyPbl 

[1] V.V. Sokolov and A.V. Turbiner, Quasi-exact-solvability of the A 2 Elliptic model: algebraic form, sl(3) hidden 
algebra, polynomial eigenfunctions, accepted to Journal Physics A, nlin. arXiv: 1409.7439 

[2] W. Riihl and A.V.Turbiner, Exact solvability of the Calogero and Sutherland models, Mod. Phys. Lett. AlO, 
2213-2222 (1995) 

[3] A.G. Matushko, On polynomial form for the elliptic Calogero-Moser system, to be published 

[4] M.A. Olshanetsky and A.M. Perelomov, Quantum integrable systems related to Lie algebras, Phys. Repts. 
94 (1983) 313-393 

[5] G. Darboux, Sur une fequation linfeare. C. R. Acad. Sci. Paris, t. XGIV (1882), no. 25, 1645-1648. 

[6] K. Takemura, Quasi-exact solvability of Inozemtsev models J. Phys. A: Math. Gen. 35 (2002), 8867-8881. 

[7] A.V. Turbiner, Quasi-Exactly-Solvable Problems and the SL(2,R) Group, Gomm.Math.Phys. 118 (1988), 
467-474 

[8] B.A. Dubrovin, Geometry of 2D topological field theories. In Integrable Systems and Quantum Groups, 
Lecture Notes in Math. 1620 (1996), 120-348. 

[9] A.V. Turbiner, Lie-algebras and Linear Operators with Invariant Subspaces, in Lie Algebras, Cohomologies 
and New Findings in Quantum Mechanics, AMS ‘Contemporary Mathematics’ 160 (1994), 263-310. 

[10] A. P. Veselov On Darboux-Treibich-Verdier potentials, Lett. Math. Phys. 96 (2011), no. 1-3, 209-216. 

[11] V.l. Inozemtsev, Lax representation with spectral parameter on a torus for integrable particle systems, Lett. 
Math. Phys. 17 (1989), 11-17. 


